Constraining Hq from Lyman-a Forest and Baryon Acoustic 

Oscillations 



V. C. Busti 

Departamento de Astronomia, Universidade de Sao Paulo, USP, Sao Paulo, SP, Brazil 

05508-090 

vcbusti@astro . iag.usp.br 

R. N. Guimaraes 

Programa de Modelagem Computacional, SENAI - Cimatec, Salvador, BA, Brazil 

41650-010 

rguimara@eso . org 

and 
J. A. S. Lima 

Departamento de Astronomia, Universidade de Sao Paulo, USP, Sao Paulo, SP, Brazil 

05508-090 

limaj as@astro . iag . usp . br 



Received 



accepted 



- 2 - 
ABSTRACT 

A new method is proposed to measure the Hubble constant H through the 
mean transmitted flux observed from high redshift quasars. A semi-analytical 
model for the cosmological-independent volume density distribution function is 
adopted which allows one to obtain constraints over the cosmological parameters 
once a moderate knowlegde of the InterGalactic Medium (IGM) parameters is 
assumed. By assuming a flat ACDM cosmology we show that such method alone 
cannot provide good constraints on the pair of free parameters (h, fl m )- However, 
it is possible possible to break the degeneracy on the mass density parameter by 
applying a joint analysis involving the baryon acoustic oscillations (BAOs). Our 
analysis based on two different samples of Lyman-a forest resctricts (h, fl m ) on 
the intervals 0.58 < h < 0.91 and 0.215 < tt m < 0.245 (la). Although the 
constraints are weaker comparatively to other estimates, we point out that with 
a bigger sample and a better knowledge of the IGM this method may present 
competitive results to measure the Hubble constant independently of the cosmic 
distance ladder. 

Subject headings: cosmological parameters — distance scale — intergalactic medium 
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1. Introduction 

Measurements of the Hubble constant Hq are vital to the establishment of a cosmic 
concordance model in cosmology. Such constant plays a role in most of cosmic calculations, 
as the age of the Universe, its size and energy density, primordial nucleosynthesis, physical 
distances to objects and so on. Its importance must even increase in the next decade 
because new missions and observational projects are being designed to provide accurate 
measurements of the Hubble constant thereby also increasing the precision in the other 
cosmological parameters. 

There are several ways to determine H Q by using low and high redshift sources. The 
most commom methods are based on the period-luminosity relation seen in Cepheids, tip of 
the red giant branch and Supernovae observed in the local Universe (Freedman & Madore 
2010). However, some alternative procedures like the Sunyaev-Zel'dovich effect combined 
with X-ray emission from clusters allow a measurement of H from high-redshift objects 
have also been discussed (Cunha et al. 2007, Holanda et al. 2010). Although considering 
that some errors recently reported in the literature are small (e.g. Riess et al. 2009), 



Sandage et al. 



(120061 ) 



discrepancies among different measurements as the value given by 
show that independent estimates are needed in order to achieve a reliable value for H not 
plagued by systematic errors from astrophysical environments, as well as from calibrations 
associated to the cosmic distance ladder. In this concern, the importance and cosmological 
interest on different estimates of H Q independent of the distance ladder has also been 
discussed by many authors (see Jackson 2007 and Refs. therein for a recent review). 

A possible alternative procedure to measure H is to consider the Lyman-a (Lya ) 
forest data that probes the low density intergalatic medium (IGM) over a unique range of 
redshifts and environments as seen in the spectra from high redshift quasars. The observed 
mean transmitted flux depends on the local optical depth, which in turn depends on the 
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expansion rate of the Universe H(z) and some properties characterizing the intervening 
medium along our line of sight to any quasar. In principle, with a moderate knowledge 
of the IGM properties, as the hydrogen photoionization rate and mean temperature, the 
values of some cosmological parameters including the Hubble constant can be constrained. 

In the last decade, several studies based on Lya forests as cosmological tools have been 
performed (Macdonald et al. 2000, Weinberg et al. 2003, Wyithe et al. 2008, Viel et al. 
2009). Some works use the flux power spectrum of the Lya forest to infer the cosmological 



parameters in two ways: inversion o 



matter power spectrum rtCroft et al 



use the power spectrum directly (IMcDonald et al. 



the f lux power 



1998 



1999 



H 



ui 



spect r um to get the underlying dark 



1999; 



2000 



Nusser &; Haehnelt 



2005 



Mandelbaum et al. 



1999), or to 



2003 ). 



On the other hand, a simpler approach is t o ado pt a semi-analytical model to describe 



the IGM as proposed by 



Miralda-Escude et al 



(120001 ). who derived a fitting formula for the 



volume density distribution function in agreement with simulations which provided si milar 



19971 ). 



results for different numerical methods and different cosmologies (IRauch et al. 
Therefore, a theoretical mean transmitted flux is obtained and compared to observational 
data in order to constrain the cosmological parameters. 

In this letter, by assuming a flat ACDM cosmology, we apply the above described 
procedure for deriving new constraints on the Hubble constant H and the matter density 
parameter Q m . By considering two independent samples of Lya forest as compiled by 
Bergeron et al. (2004) and Guimaraes et al. (2007), respectively, two distinct statistical 
analyses are performed. Firstly, we consider only the Lya forest data set whereas the 
second approach involves a joint analysis combining the Lya forest data with the SDSS 



measurements o 



( jEisenstein 



the baryon acoustic peak from the clustering of luminous red galaxies 



20051 ). As we shall see, the degeneracy on the free parameters defining the 



(h, Q m ) plane is naturally broken when the Lya data are combined with the standard ruler 
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as given by the baryon acoustic oscilations thereby providing a new independent method to 
estimate the Hubble parameter. 



2. The Lya Forest Samples 



As remarked earlier, in the present analysis we consider two different samples. The 
first one consists of 18 high-resolution high signal-to-noise ratio (S/N) spectra. The spectra 
were obtained with the Ultra- Violet and Visible Echelle Spectrograph (UVES) mounted on 
the ESO KUEYEN 8.2 m telescope at the Paranal observatory for the ESO-VLT Large 
Programme (LP) 'Cosmological evolution of the Inter Galactic Medium' (PI Jacqueline 
Bergeron; Bergeron et al. 2004). The spectra were taken from the European Southern 
Observatory archive and are publicly available. This programme has been devised to gather 
a homogeneous sample of echelle spectra of 18 QSOs, with uniform spectral coverage, 
resolution and signal-to-noise ratio suitable for studying the IGM. The spectra have a 
signal-to- noise ratio of 40 to 80 per pixel and a spectral resolution A/ A A ~ 45000 in 
the Lya fore s t regi on. Details of the proce dure used for data reduction can be found in 



Chand et al 



( 120041 ) and 



Aracil et al. 



( 120041 ) 



The second sample consists of 45 qua sars obtained 



the spectra were measured with the ESI (ISheinis et al. 



Guimaraes et al. 



(120071 ). where 



20021 ) mounted on the Keck II 



10m telescope. The spectral resolution is R ~ 4300 with a signal-to-noise ratio usually 
larger than 15 per 10 km s _1 pixel. This sample has also been studied Prochaska and 
collaborators (2003a,b) in their discussions related to the existence of Damped Lya systems. 
In the present analysis we have used the Lya forest over the rest-wavelenght range 
1070-1170 A to avoid contamination by the proximity effect close to the Lya emission line 
and possible OVI associated absorbers. We also carefully avoided regions flagged because 
of damped/sub-damped Lya absorption lines. We divided each spectrum in bins of length 
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100 A in the observed frame, corresponding to about Az = 0.08, for the last bin we have 
used a lenght of 500 A to keep the statistical homogeneity. The mean transmitted flux was 
calculated as the mean flux over all pixels in a bin. At each redshift, we then averaged the 
transmitted fluxes over all spectra covering this redshift. 

In Fig [TJ we display the mean transmitted flux for the two samples. Notice that all 
errors were estimated as the standard deviation of the mean values divided by the square 
root of the number of spectra. The corresponding scatter in the transmitted flux cannot 
be explained by photon noise or by uncertainties in the continuum thereby suggesting that 
errors are dominated by cosmic variance. 



3. Basic Equations 

The basic measurement obtained from the Lya forest is the mean transmitted flux at 
redshift z defined as 



(F)(z) = (e W (-r(z))), 



(1) 



where r(z) is the local optical depth a nd angle brackets de note an average over the 



line of sight. Following standard lines (IHui fc Gnedin 



19971 ). throughout this paper 
photoionization equilibrium and a power-law temperature-density relation, T = T (l + 8)?, 
for the low-density IGM will be assumed, where 8 stands for the local overdensity and T 



for the 



( [Peebles 



GM temperature at mean density (8 = 0). The 



1993 



Padmanabhan 



2002 



Faucher-Giguere et al. 



ocal op tical depth is given by 



2008aJ) 



■(z)=A(z)(l + 8) 



2-0.75,9 



(2) 



with 
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A(z) = 



ve 2 f Lya f PcTit Q b V X(X + 0.5Y) /,'„•/„ 



1 + ^) 



m e v Lya \ m p j H(z) V 

In the expression above, is the oscillator strength of the Lya transition, v^ yOL 
is its frequency, m e and m p are the electron and proton masses respectively, fib is the 
baryon density parameter, H(z) is the Hubble parameter, p cr - lt is the critical density, X 
and Y are the mass fra ctions of hydrogen and helium respectively taken to be 0.75 and 



(3) 



0.25 (Buries et al. 



20011 ). #0 = 4.2 x 10~ 13 cm 3 s-V(10 4 K)~ a75 and V is the hydrogen 



photoionization rate. 

The mean transmitted flux is obtained by integrating the local optical depth through a 
volume density distribution function for the gas A = 1 + 8, 



(F)(z) 



POO 

/ dAP(A;z)exp[-r(z)]. 
Jo 



(4) 



As is widely known, 



Miralda-Escude et al. 



(120001 ) derived an approximate analytical 



functional form for the distribution function given by 
where the parameters A and Cn are derived by re quiring the total volume and mass to be 



(5) 



normalized to unity. 



Miralda-Escude et al. 



(120001 ) extrapolated the distribution function 



and obtained 5q = 7.61/(1 + z) wit h an accuracy better t han 1 % from fits to a numerical 



simulation at z =2, 3, and 4 from 



Miralda-Escude et al 



(119961 ). In order to apply this 



formalism to our data, the values of b were derived from a cubic interpolation from the 
values in the simulation. Applications of this distribution function to constrain cosmological 
parameters is well justified since simulatio ns with different n umerical methods and different 



cosmologies provided very similar results ( jRauch et al. 



19971 ). 



From now on a flat ACDM cosmology will be assumed. In this case, the Hubble 
parameter is given by 
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h(z) = h [n m (i + z f + (i - n m )] 1/2 , (6) 

where we will adopt the convention h = H /100 km s -1 Mpc -1 , which is the Hubble 
constant normalized in units of 100 km s _1 Mpc -1 . 

The effects of H in the mean transmitted flux are sho wn in Fig. [TJ along with our 
observational data (red circles) and 



Guimaraes et al. 



(120071 ) (black squares). We see that 



with our sample we ob tain a scattered data set comparatively to the sample used by 



Guimaraes et al 



(120071 ) . which is related to the sample size. It is also evident that our data 
cover most of values for the h parameter, so we do not expect tight constraints from the 
Lya data set alone. 



4. Analyses and Results 

Let us now perform a statistical analysis to find the constraints on the cosmological 
parameters. The full set of parameters are represented by p = (h, Q m , T , T n , flhh 2 , 13). We fi x 



the value of f2 fc /i 2 = 0.0218 using the latest observations of de uterium (jPettini et al 



from the Big Bang Nucleosynthesis rtSimha &: S tcigman 



is also considered which gives (3 = 0.62 ( )Hui fc Gnedin 



20081 ) 



2008J). An early reionization model 



19971 ). The posterior probability of 



the parameters P(p\d) given the data d is 



P(p\d) 



P(d\p)P(p) 
P{d) 



(7) 



where P{d) is a normalization constant, P(p) is the prior over the parameters and 
P(d\p) oc e-x 2 /2 is 

the likelihood with the usual definition, for the Lyo; forest data, of 
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^Lya 



((F lh )(z i; p)-{F)(z,)y- 



a 



(8) 



where (i^)(^;p) is the theoretical mean transmitted flux, (F)(zi) is the observational 
mean transmitted flux and <J(F)(zi) is its respective uncertainty. We treat the combination 
TTq ' 75 as a nuisance parameter with a flat p rior. The ranges chosen are T = [2, 2.5] x 10 4 



K from estimates of 



Zaldarriaga et al. 



measurements reportec 



2001 



Meiksin fc White 



in the literature ( 



2004 



fl200lh and T =J0-8.1 



?lauch et al. 



Tytler et al 



2004 



1997 



x 10 s which cover some 



McDonald & Miralda-Escude 



Bolton et al 



although in disagreement with the values obtained by 



2005 



Kirkman et al 



Faucher-Giguere et al 



2005), 



fl2008a 



b) 



who used a redshift-dependent relation for T not favoured by our data. 

In what follows, we first consider the Lya forest data separately, and, further, we 
present a joint analysis including the BAO signature extracted from the SDSS catalog 
(Eisenstein et al. 2005). 



4.1. Limits from Lyct Forest data set 

In Fig. [2] we display the results of the statistical analysis performed with the Lya forest 
data. We see that the data do not provide good constraints to both parameters. Several 
sanity checks were performed in order to determine the values fixed in the statistical 
analysis. In particular, we have seen that different values for /3, different intervals for 
py- •< j n the marginalization, or even a redshift-dependence for To, have resulted in 
poorer fits compared to what is shown in figure [2j In addition, even considering that 
cosmological constraints are affected by the IGM physical parameters the opposite vision is 
also true, that is, different cosmologies may also provide results not compatible for the IGM 
parameters. In this context, it is recommended to fit cosmological and IGM parameters 
together in order to obtain more trustful results. 
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Despite the considerations discussed above, we see from Fig. [2] that the slope in the 
(h, Q m ) plane suggests that a joint analysis with an independent test constraining only fl 7 
could provide interesting limits to the Hubble parameter. Thus, a joint analysis with the 
BAO data is presented in the next subsection. 



4.2. Lya Forest and BAO: A Joint Analysis 



In order to achieve better constraints on the cosmic parameters we apply a joint 
analysis involving Lja forest data and baryon acoustic oscillations data obtained from 



46748 luminous red galaxies se. 



represented by the parameter (jEisenstein 



ected from the SPSS Main Sample. The BAO scale can be 



2005|) 



A = 



n 1/2 



H{z* 



,1/3 



-T(z ± 



2/3 



0.469 ±0.017, 



(9) 



where = 0.35 is the redshift at which the acoustic scale has been measured, 'H(z) is 
H(z)/H and F(z*) is the dimensionless comoving distance to z*. 

From equation ([9]) it is seen that the BAO scale is independent of h thereby yielding 
constraints only for the matter density parameter. The statistical analysis is performed 
with x 2 = xlya + Xbao, where 



Xbao 



A - 0.469 
0.017 



(10) 



In Fig. [31 we show the limits on the pair of parameters (h, Q m ) obtained from our 
joint analysis involving Lya forest data and BAO. Within 68.3% confidence level, we have 
obtained 0.215 < Q m < 0.245 and 0.58 < h < 0.91. The value for Q m is lower than what 
was found by other observations (e.g. 



Amanullah et al. 



20101 ). but it is likely that the 
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difference is due to systematic errors not taken into account in our analysis. 



In Tab. (TJ we show some recent measurements of H using di fferent techniques an d 



data. Although existing a tension among different estimates (e.g. 



Paraficz fc Hjorth 



Sandage et al. 



2006 



20101 ) . our approach do not provide stringent limits which would allow a 



decision about the correct value of H . The interesting aspect is that the method discussed 
here provides an independent estimate which is in agreement with the values obtained 
adopting different approaches (see Tabled]). 



5. Conclusions 

In this work we have used a cosmological-independent semi-analytical model to 
describe the IGM and data from Lya forest and baryon acoustic oscillations to constrain 
cosmological parameters. Limits on the Hubble constant Ho and the matter density 
parameter Q m were derived by assuming a flat ACDM cosmology. 

By applying a statistical analysis using only the Lya forest data we did not obtain 
good constraints to both parameters. However, when we applied a joint analysis involving 
the Lya forest data and BAO interesting limits were found, being the parameters restricted 
to the intervals 0.215 < Vt m < 0.245 and 0.58 < h < 0.91 within 68.3% confidence level. 
These results are in agreement with recent measurements reported in the literature (see 
Tab. [T]), but are weaker due to the limited sample and our poor knowledge of the IGM. 
We expect better constraints to the Hubble constant in the near future by using a bigger 
sample and a deeper understanding of the IGM, which will turn this technique competitive 
with other estimates with the advantage of being independent of a cosmic distance ladder. 

The authors are grateful to J. V. Cunha and J. F. Jesus for helpful discussions. VCB is 
supported by CNPq and RNG is partially supported FAPESB. JASL is partially supported 
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Fig. 1. — Mean transmitted flux as a function of redshift for Vt m = 0.23, T = 0.8 x 10 s 
T = 2.5 x 10 4 K and some selected values of the h parameter. The data points correspond 



to the Lya forest measurements obtained by this work (red circles) and 
(120071 ) (black squares). 



Guimaraes et al. 
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Fig. 2. — Confidence regions (68.3%, 95.4% and 99.7%) in the (fl m , h) plane provided by 
Lya forest data. The basic conclusion is that Lya forest alone cannot constrain the pair of 
parameters. 
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Fig. 3. — Contours in the Q m — h plane provided by the joint analysis combinig Lya Forest 
and BAO. As before, the contours correspond to 68.3%, 95.4% and 99.7% confidence levels. 
Note that the best-fit model converges to h = 0.67 and fl m = 0.23. These results should be 
compared with the ones presented in Figure 2. 



-20 - 



Table 1: Limits to H 



Method 



Reference 



Cepheid Variables 

Age Redshift 

SNe Ia/Cepheid 

SZE+BAO 

Old Galaxies + BAO 

SNe Ia/Cepheid 

CMB 

SZE+BAO 
Time-delay lenses 
Lyo + BAO 



Freedman (2001 ) (HST Project) 
Jimenez et al. (2003) (SDSS) 
Sandage et al. (2006) 
Cunha et al. (2007) 
Lima et al. (2009) 
Riess et al. (2009) 
Komatsu et al. (2011) (WMAP) 
Holanda et al. (2010) 
Paraficz fc Hiorth (2010) 
This work 



0.72 ± 0.08 

0.69 ±0.12 

0.62 ± 0.013(rand.)±0.05(syst. 
n 74+0.04 

u - '^-0.03 

0.71 ± 0.04 
0.742 ± 0.036 
0.710 ±0.025 



n 7QO+0.043 
u - '°^-0.037 



0.76 ±0.03 



0.67 



-0.24 
-0.09 



